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Metal-insulator transitions in GdTiOa/SrTiOs superlattices 
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The density functional plus U method is used to obtain the electronic structure, lattice relaxation 
and metal-insulator phase diagram of super lattices consisting of m layers of Gadolinium Titanate 
(GdTiOs) alternating with n layers of Strontium Titanate (SrTiOa). Metallic phases occur when 
the number of SrTiOa layers is large or the interaction U is small. In metallic phases, the mobile 
electrons are found in the SrTiOa layers, with near-interface electrons occupying xy-devWed bands, 
while away from the interface the majority of electrons reside in xzfyz bands. As the thickness 
of the SrTiOa layers decreases or the on-site interaction U increases a metal-insulator transition 
occurs. Two different insulating states are found. When the number of SrTiOa layers is larger than 
one, we find an insulating state with two sublattice charge and orbital ordering and associated Ti-0 
bond length disproportionations. When the number of SrTiOa units per layer is one, a different 
insulating phase occurs driven by orbital ordering within the quasi one-dimensional xzfyz bonding 
bands connecting Ti atoms across the SrO layer. In this phase there is no sublattice charge ordering 
or bond disproportionation. The critical U for the single-layer insulator is ~ 2.5 eV, much less than 
critical U ~ 3.5 eV required to drive the metal-insulator transition when the number of SrTiOa is 
larger than one. Inconsistencies between the calculation and the experiment suggest that many-body 
correlations may be important. A local inversion symmetry breaking around Ti atoms suggests the 
possibility of in-plane ferroelectric polarization in the insulating phase. 

PACS numbers: 73.20.-r, 71.30.+h, 75.70.-i, 71.27.+a 


I. INTRODUCTION 

Complex oxide systems have been investigated exten¬ 
sively because of the variety of interesting phases that 
may be achieved by alloying and pressure^. Of partic¬ 
ular interest are materials with structures based on the 
ABO3 perovskite motif (here A represents a rare earth 
or alkali atom, B a transition metal, and O an oxygen), 
because the many available possibilities for the A and B 
site ion lead to a great diversity of interesting properties. 
Over the last decade it has become clear that complex ox¬ 
ide heterostructures may be grown with precise layer-by- 
layer control of the composition^“^ enabling the creation 
of what are in effect new materials with potentially new 
properties. Moreover, carrier doping by charge transfer 
across interfaces can be much higher than that obtainable 
by alloying or vacancy formation in bulk materials^’^’^ 
while electronic gating allows the continuous control 
of carrier doping.Heterostructures comprised of 
two different ABO3 transition metal oxide perovskites 
with comparable bulk lattice constants have been re¬ 
ported to exhibit metal-insulator transitions (MIT)^^^^^, 
magnetism^^ as well as superconductivity.^’^^ 

The electronic properties of semiconductor het- 
erostructures are determined by bulk band gaps, work 
functions, and the position of donor (acceptor) levels. 
In many transition metal oxides, electronic properties 
are, in addition, sensitive to structural features including 
octahedral rotation angles and transition-metal-oxygen 
bond lengths. Thus local structural changes across in¬ 
terfaces may play an important role in determining the 
electronic properties of oxide heterostructures. Moreover 
structural distortions propagate only a few lattice con¬ 
stants, so superlattice thickness has a structural as well 


as a quantum confinement effect. Undertanding the in¬ 
terplay of the various factors that contribute to metal- 
insulator transitions in oxide superlattices is an impor¬ 
tant open problem. 

Heterostructures comprised of Mott-insulating 
GdTi03(GT0) and band insulating SrTi03(ST0) 
provide an interesting model system. A single interface 
separating semi-infinite slabs of GdTi03 and SrTi03 is 
found to be metallic with sheet charge density about 
a half electron per in-plane unit cell,^ consistent with 
elementary ‘polar catastrophe’ notions of interface 
doping and in contrast with the more complicated 
behavior of the widely studied LaA103(LA0) systems, 
where the density of mobile carriers is much less than 
the polar catastrophe value^^. Part of the issue may 
be that GdTi03 and SrTi03 share a common Ti—O 
network so that the main defects at an interface would 
be Gd—Sr antisite defects, which might be expected 
to disrupt the electron structure less than the Al—Ti 
antisite defects that might occur at the LAO-STO 
interface. A further interesting feature is that GdTi03 
is itself a Mott insulator, and has a large amplitude 
rotational distortion away from the basic cubic ABO3 
perovskite structure. Taken together, these features 
suggest that the GTO—STO system may exhibit an 
interesting interplay of structural, interface, quantum 
confinement, and correlation effects, perhaps not too 
badly complicated by disorder, making it a suitable 
model system for investigation of general issues relating 
to metal-insulator transitions in oxide heterostructures. 

Recently, thickness-dependent metal insulator 
transitions have been reported in GTO—STO 
heterostructures.^ As the thickness of SrTi03 de¬ 
creases to two unit cells of SrTi03, the interface is 
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found to become insulating. Both in metallic and in 
insulating interfaces, variations of the crystal structure 
that decay within several unit-cells of the interface are 
observed.In metallic interfaces, quantum oscillation 
measurements indicate that the conducting carriers 
move primarily in the plane of the interface.The 
theory of the insulating phase was investigated on the 
basis of first-principles DFT+U calculations by Chen 
and Balents,^^ who proposed a novel dimerization 
mechanism that could lead to an insulating ground state 
when the number of SrTiOs layers, n = 1 and by Bjaalie 
et al?^ who reported an insulating ground state with 
n = 2 using DFT with a hybrid functional. In this 
paper we present a more general DFT+U investigation 
of the metal-insulator phase diagram and structural 
properties of the GTO-STO system. Although DFT+U 
is a Hartree approximation which does not capture the 
full complexity of many-body physics, it enables the 
investigation of the interplay between structural relax¬ 
ation and electron correlation effects. The changes in 
the electronic structure in different structural phases are 
investigated and the spatial distribution of the electron 
gas and its orbital polarization for different thickness 
of SrTiOs are determined as function of interaction 
strength, U. 



FIG. 1: (color online) Schematic phase diagram of electronic 
phases as a function of SrTiOs layer thickness n and on-site 
Coulomb repulsion U. 


Our crucial finding, illustrated in Fig. 1, is that there 
are two different insulating phases, one occurring for 
ultra-thin (n = 1) and one for thicker (n > 1) STO lay¬ 
ers. In the DFT+U approximation employed here, some 
form of charge, orbital, or magnetic ordering is required 
to stabilize insulating phases. When there is more than 
one SrTiOs unit cell (n > 1), we find a MIT from ferro¬ 
magnetic metal to charge-ordered insulator with a nar¬ 
row intermediate regime of charge ordered low density of 
states behavior, which we interpret as the signature in 
this calculations of a “bad metal” phase.The metal- 
insulator transition is accompanied by a structural transi¬ 
tion resulting in sublattice bond disproportionation, and 


within DFT+U the charge ordering is the main driving 
force for insulating/bad-metallic behavior. In the case 
of a single SrO layer (n = 1) we find two insulating 
phases. In one of them, yzjxz orbitals on opposite sides 
of the SrO layer form bonding and antibonding states; 
the antibonding states are empty and the spin-polarized 
half-filled bonding bands become gapped by orbital or¬ 
dering. This transition from ferromagnetic metal to fer¬ 
romagnetic insulator phase is consistent with the pre¬ 
vious Hartree-Fock calculations of Chen and Balents.^^ 
The critical interaction strength required to drive this in¬ 
sulating phase is significantly lower than that needed to 
drive the insulating phase I. As on-site Coulomb interac¬ 
tion increases at n = 1 a transition occurs to a different 
insulating phase, characterized by two-sublattice charge 
order and similar to the insulating phase found at n > 1. 

The rest of the paper is organized as follows. In sec¬ 
tion II we state the calculation methods that we use. 
Section III presents structural phase transitions in the 
different crystal structures. Electronic structures in each 
phase are are presented in Section IV. Section V is de¬ 
voted to an analysis of the electronic structure with a 
single SrO layer. In Section VI we present structural and 
electronic phase diagram in terms of layer thickness and 
electron correlation. Section VII presents possible fer¬ 
roelectric polarization from broken inversion symmetry 
and centrosymmetry. In Section VIII we discuss a sim¬ 
plified view of metal insulator transition based on band 
alignment, followed by summary in Sec, IX. 


II. CALCULATIONAL METHODS 

The electronic structure and atomic structures are cal¬ 
culated using the GGA-fU method as implemented in the 
Vienna ab-initio simulation package (VASP).^^’^^ We use 
a plane wave basis set with energy cutoff 500 eV and the 
projected augmented wave method. For all superlattices 
we use a unit cell \f2 x \f2 x 2/ (/ : integer) structure to 
describe proper octahedral rotations and orbital ordering 
of the superlattices. All the structures are fully relaxed 
while the in-plane lattice constant is constrained to the 
substrate {a = b = 3.86A). Gonvergence is reached if 
the energy difference between two consecutive iterations 
is within 0.1 meV for electronic iterations and 1 meV 
for ionic relaxations. The minimum /c-point grid size is 
6x6x2. For on-site Coulomb interaction the rotation- 
ally invariant form^^ is used for Ti-d orbitals. We treat 
the Hubbard U as an adjustable parameter while fixing 
the value of J as 0.6 eV. 


III. LATTICE STRUCTURES 

In this section we investigate the lattice structures of 
(GTO)m(STO)n superlattices for different n and differ¬ 
ent electronic phases. We begin with a brief discussion 
of the bulk structures. For SrTiOs we focus on the high 
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FIG. 2: (color online) Structural properties of (GTO) 4 (STO )4 
superlattices in the metallic phase of the superlattice com¬ 
puted as described in the main text with f/ = 3 eV. (a) The 
atomic structure of the super lattice. A and B represent two 
inequivalent in-plane Ti atoms and the numbers in the bot¬ 
tom label the Ti02 planes, (b) Octahedral volumes around 
A and B sublattice Ti atoms. The green dashed and blue 
dotted lines are respectively the octahedral volumes of bulk 
GTO and STO calculated with experimental in-plane lattice 
constant, (c) Average values of in-plane (^||) and out-of-plane 
{0±) Ti-O-Ti angles. The green dashed and blue dotted lines 
are the in-plane and out-of-plane octahedral angle calculated 
for bulk GTO with experimental substrate, respectively, (d) 
In-plane (x, y) Ti-0 bond lengths of A sublattice Ti atoms, 
(e) Out of plane (z) bond lengths of A sublattice Ti atoms. 
The six oxygen atoms around the each Ti atom are labeled by 
the direction from Ti to oxygen atoms expressed by x±,y±, 
and where the direction {x,y, z) and subscripted sign rep¬ 
resent pseudo-cubic axis dehned in panel (a) and plus and 
minus direction, respectively. Exchanging x± to y^p gives the 
B sublattice Ti-0 distances. 


temperature phase which is nearly cubic with a crystal¬ 
lographic unit cell containing one Ti ion and an octa¬ 
hedral volume of ^ 9.9 A^. (In bulk STO, a transition 
occurs at T ^ IlOK to a lower symmetry phase only 
slightly distorted from the high temperature structure. 
The differences between the high and low T phase are 
not important for our discussion.) GdTiOs is strongly 
distorted away from the cubic perovskite structure, with 
a unit cell containing two Ti ions, a Ti-O-Ti bond an¬ 
gle ^ 145° and a substantially increased octahedral vol¬ 
ume ^ 11.3 A^. Further GdTiOs displays a high degree 
of (7r,7r,0) “orbital ordering”, involving a two sublattice 
pattern in which in a plane (which we take to be perpen¬ 
dicular to the growth direction) the Ti-0 bonds in the 
X direction are elongated on one Ti sublattice and com¬ 
pressed on the other, with the bonds in the y direction 
behaving oppositely. The differences in the octahedral 
rotations will be seen to have important consequences. 
To capture these we choose a \/2 x \/2 x 2/ (/ : integer) 
computational unit cell such that the two inequivalent Ti 
ions lie in the plane perpendicular to the growth direc¬ 
tion. 

Fig. 2 presents relaxed lattice structures obtained using 
a moderate U = 3 eV so that the superlattice is metallic. 
The top panel (a) gives atomic positions. The second 
panel (b) shows that the octahedral volume changes dis- 
continuously across the interface, taking essentially the 
bulk GdTiOs value for Ti ions bounded on all sides by Gd 
ions and taking essentially the bulk SrTiOs value for Ti 
ions either bounded on all sides by Sr ions, or at the in¬ 
terfaces with two Gd and two Sr neighbors. The spatial 
variation of the structure is related to the band align¬ 
ments, which are such that the charge density for all Ti 
ions with eight Gd neighbors takes essentially the bulk 
GdTiOs value. We will show below that the correlation 
effects expressed here by the U make an important con¬ 
tribution to the band alignments. The third panel (c) 
shows the variation of the Ti-O-Ti bond angle for both 
bonds in the planes (6>||) and along the growth direction 
{0±). We see again that the bond angle deviates sig¬ 
nificantly from the bulk GdTiOs value only in Ti layers 
surrounded by Sr or for “apical” oxygens connecting an 
interface Ti to a Sr. 

The remaining panels show Ti-0 bond lengths. Bulk 
GdTiOs has a two-sublattice orbital ordering; our com¬ 
putational unit cell is such that the ordering wavevec- 
tor is perpendicular to the growth direction (as is phys¬ 
ically reasonable), and the resulting order is an in-plane 
xzjyz alternation. Fig. 2 (d) shows the in-plane bond 
lengths for one of the two sublattices. We see that in 
the GdTiOs region (layers 1,2,8) the Ti-0 bond length 
in the plus and minus x directions is ^ 0.1 A greater 
than the Ti-0 bond length in the plus and minus y di¬ 
rections (on the other sublattice, the pattern is reversed). 
Similarly within the SrTiOs (layers 4,5,6) there is neg¬ 
ligible bond-length disproportionation. Thus as far as 
the in-plane bonds are concerned, the lattice distortions 
associated with GdTiOs propagate only as fas as the in- 
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FIG. 3: (color online) Structural properties of (GTO) 4 (STO )4 
superlattices in the insulating phase computed as described 
in the main text with = 4 eV. Atomic coordinates are 
same as shown in Fig. 2. (a) Octahedral volumes around 

A and B sublattice Ti atoms. The green dashed and blue 
dotted lines are respectively the octahedral volume of bulk 
GTO and STO calculated with experimental in-plane lattice 
constant, (b) Average values of in-plane (^||) and out-of-plane 
Ti-O-Ti angles. The green dashed and blue dotted lines 
are the in-plane and out-of-plane octahedral angle for GTO 
calculated with experimental substrate, respectively, (c-f) Ti- 
O distance between A and B sublattice Ti atoms following the 
same definition in the Fig. 2. 


ter face layer, consistent with what was observed for the 
bond angles and octahedral volumes. 

However, the in-plane bonds in the interface layer and 
the out-of-plane bonds throughout the SrTiOs region 
show an interesting behavior implying breaking of the 
local inversion symmetry at the Ti site. The inversion 
symmetry breaking in the 2 : direction may be under¬ 
stood in terms of the dielectric properties of SrTiOs. 
The electric fields associated with the polar discontinuity 
of the GdTiOs/SrTiOs interface and the associated in¬ 
duced charges are partly screened by polar modes of the 
SrTiOs lattice (in particular an off-centering of the Ti) 
as discussed in Okamoto et . Our new finding is a 
difference in the x and y bonds implying a spontaneous 
ferroelectric distortion in the plane. We also note that 
there is no difference in the octahedral volumes between 
A and B sublattices implying equal electron occupancy of 
the A and B sublattice Ti states. We define this structure 
with the same octahedral volume for A and B sublattices 
as the non charge-ordered (NCO) phase. 

In Fig. 3 we present the structural properties of insulat¬ 
ing (GT0)4(ST0)4 superlattices obtained using DFT+U 
with U = A eV. The variation of octahedral volume is 
very similar to that found in the previous metallic case, 
except that in the interface layer a clear difference be¬ 
tween sublattices is visible. One sublattice (“A”) has 
almost the same octahedral volume as in bulk GdTiOs 
while the other has almost the same octahedral vol¬ 
ume as bulk SrTiOs. This difference reflects the almost 
complete charge order noted by many previous stud¬ 
ies of related interfaces including LaAlOs/SrTiOs^^ and 
LaTiOs/SrTiOs-^^ The difference in octahedral volumes 
between sublattices has almost no effect on the in-plane 
bond angles; however the out of plane bond angels involv¬ 
ing apical oxygens around interface Ti ions show small 
differences (about 3°) between A- and B- sublattices 
(not shown here.) The magnitude of the centrosymmetry 
breaking distortions in the ^-direction Ti-0 bond is much 
less, reflecting the almost complete confinement of polar 
charge to the interface layer, strongly reducing the need 
to screen internal electric fields. We define this structure 
with inequivalent octahedral volumes as the charge order 
(CO) phase. 

For general (GTO)m(STO)n superlattices, the struc¬ 
tural transition from the NCO to the CO phase is driven 
by increasing the on-site Coulomb interaction U. How¬ 
ever for n = 1 superlattices with a single STO unit cell, 
the distortion around the interface is different when U is 
small. Fig. 4 shows the structure for the metallic inter¬ 
face found inam = 5,n = l superlattice with U = 2 
eV. As in the n > 1 case there is no difference in the 
octahedral volume between A and B sublattice but there 
is a small difference in octahedral volume between layer 
3 and 4 mainly due to the difference in the Ti-0 distance 
in 2 ; direction. This is different from the NCO phase of 
n > 1 superlattices in which the octahedral volumes of 
two interfaces are different. We define the structure as 
CO H phase. We also note that the amount of Jahn- 
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FIG. 4: (color online) Structural properties of (GT0)5(ST0)i 
superlattices computed as described in the main text with 
U = 2 eV. (a) The atomic structure of the superlattice. A 
and B represent two different in-plane Ti atoms and the num¬ 
bers in the bottom denote Ti02 plane number, (b) Octahe¬ 
dral volumes around A and B sublattice Ti atoms. The green 
dashed and blue dotted lines are respectively the octahedral 
volume of bulk GTO and STO calculated with the experi¬ 
mental lattice constant, (c) Average values of in-plane (^||) 
and out-of-plane {0±) Ti-O-Ti angles. The green dashed and 
blue dotted lines are the in-plane and out-of-plane octahe¬ 
dral angle for GTO calculated with experimental substrate, 
respectively, (d-e) Distance between Ti and O in the in-plane 
(d) and growth direction (e) pseudo cubic direction of A sub¬ 
lattice Ti atoms. Exchanging x± to gives B sublattice 
Ti-0 distances. 


Teller distortion is much less in the GTO region due to 
the decreased correlation. As U increases, the structure 
changes to sublattice charge order (CO) phase as pre¬ 
sented in Fig. 5. However unlike the n > 1 case the pre¬ 
ferred orbital character is not switched between GTO and 
the larger interface octahedron since the Ti-0 distance in 
X (y) direction are longer than y (x) for all octahedra in 
A (B) sublattice. 

Figure 6 presents the critical electronic correlation Uc 
required to derive the structural phase transition between 
two distinct distortion patterns for (GTO)^(STO)n su¬ 
perlattices. Uc is defined as the value of U in which the 
total energy with CO phase becomes lower than that with 
NCO (CO II) phase. The Uc is relatively insensitive to 
the thickness of GTO and we can see a general trend that 
Uc increases as the thickness of STO increases for n > 1. 
However, there is drastic decrease in the Uc for n = 1 
superlattices. Thus we can identify two different groups 
of Uc] one around 3.5 eV for n > 1 and around 2.5 eV 
for n = 1. The difference in the critical U and structural 
transitions for n = 1 and n > 1 superlattices implies 
that the change in the underlying electronic structures 
between two groups of superlattices will be different. 


IV. ELECTRONIC STRUCTURE OF n > 1 
GTO/STO SUPERLATTICES 

This section presents the electronic structures of su¬ 
perlattices with two and more SrO layers. The CO dis¬ 
tortion that we identified in the previous section plays 
an important role in localizing the electrons. The nomi¬ 
nal valence counting implies that each GdO layer donates 
one-half electron per in-plane unit cell to each adjacent 
Ti02 layer. Thus, the interface Ti02 layer is doped by 
an half electron per in-plane Ti. Therefore within band 
theory both charge and orbital ordering are required to 
obtain insulating states. Further, as shown in panels (c) 
and (f) of Fig. 3, the apical oxygen is displaced in the z 
direction, implying the xy orbital is higher in energy than 
the xzjyz orbitals. Thus, in this case the electrons are in 
two orbitally degenerate bands and insulating behavior 
requires charge, spin, and orbital order. 

Fig. 7 shows the total density of states (DOS) of a 
(GT0)4(ST0)4 heterostructure with U = 3 (no charge 
order) and 4 eV (charge order). The Fermi energy is de¬ 
fined to be = 0. The states in between -7.5 and -2.5 eV 
are oxygen p states and those near and above the Fermi 
energy are Ti-t 2 ^ derived states. The oxygen p-states 
shifts about 0.5 eV as the values of U increases from 3 
to 4 eV while the shift in energy of the Ti-t 2 ^ states is 
about 0.2 eV. There are evident changes in DOS near 
Fermi energy. For U = A eV, the system is insulating 
with small gap less than 0.1 eV. For U = 3 eV, there are 
substantial increases in the density of states at the Fermi 
level and the system is a half-metal. As expected we 
find that the insulating ground state is associated with 
a volume difference in octahedra between interface A-B 
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FIG. 5: (color online) Structural properties of (GTO) 5 (STO)i 
superlattices in the insulating phase computed as described 
in the main text with U = 3 eV. Atomic coordinates are same 
as shown in Fig. 4. (a) Octahedral volumes around A and B 
sublattice Ti atoms. The green dashed and blue dotted lines 
are respectively the octahedral volume of GTO and STO cal¬ 
culated with the experimental lattice constant, (b) Average 
values of in-plane (^||) and out-of-plane {0±) Ti-O-Ti angles. 
The green dashed and blue dotted lines are the in-plane and 
out-of-plane octahedral angle for GTO calculated with ex¬ 
perimental substrate, respectively, (c-f) Distance between Ti 
and O in each pseudo cubic direction of A and B sublattice 
Ti atoms. 



FIG. 6: (color online) Dependence of critical interaction 
strength Uc for the structural phase transition on the num¬ 
ber n of STO layers in (GTO)m(STO)n super lattices, com¬ 
puted as described in the text for different GTO thicknesses 
m (shown with different symbols indicated in Figure legend). 
For U < Uc the phase is either not charge ordered (NGO), for 
n > 1 or weakly charge ordered (G02) for n = 1. For U > Uc 
the phase is strongly charge ordered (GO) at all n. 


sublattice, while in the metallic ground state all octahe- 
dra have equal volumes. 

In order to investigate the nature of the conduct¬ 
ing interface, we present the projected density of states 
(PDOS) of Ti-t 2 ^ derived bands for [/ = 3 eV in Fig. 8. 
For the A-sublattice Ti atoms in between GdO layers 
(panel (b) and (c) labeled as A1 and A2), the density of 
states is similar to that of bulk GTO (panel (a)) with 
negligible density of state at the Fermi level and is or- 
bitally ordered {xz on A-sublattice, yz on B-sublattice). 
The sublattice orbital ordering disappears for Ti atoms at 
the interface and in between SrO layers. At the interface 
(panel (d)), the x^-derived band is dominantly occupied 
and away from the interface the occupancy of yzjxz elec¬ 
trons is gradually increasing (panel (e-f)). These char¬ 
acteristic features of the PDOS are also found in other 
conducting superlattices with n > 1. 

In Fig. 9, we present the layer resolved occupation of 
the t 2 g orbitals of (GTO) 4 (STO )4 and (GTO) 4 (STO )6 
superlattices showing that the orbital disproportionation 
persists to relatively thick superlattices. To define the oc¬ 
cupation we integrate the PDOS of Ti t 2 g orbitals from -1 
eV to Fermi energy and normalize the result so that the 
total integrated density of states equals two electrons. 
We find that for the relatively thick STO layers stud¬ 
ied the total t 2 g occupation does not vary much in the 
STO and interface region because the screening length is 
larger then the STO thickness. We believe that for both 
(GT 0 ) 4 (ST 0)4 and (GTO) 4 (STO )6 superlattices this is 
an interface driven phenomenon similar to that found in 
LaAlOs/SrTiOs interfaces, which have a dominantly oc¬ 
cupied xy orbital in the near-interface region and yzjxz 
orbitals away from the interface^^’^^. 
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FIG. 7: (color online) Density of states of (GTO) 4 (STO )4 
superlattices, (a) DOS for U = 4.0 eV. (b) DOS for U = 3.0 
eV. (c) Gomparison of spin up density of states near the Fermi 
energy. 


Fig. 10 shows the PDOS of insulating ground state 
with U = 4 eV. In this phase we find that the electrons 
induced in the STO are tied to interface Ti atoms and 
exhibit sublattice orbital and charge ordering. As in the 
conducting interfaces, the PDOS of Ti sandwiched by 
GdO layers (panel (a) and (b)) resembles that of bulk 
GTO. At the interface, the Ti states of the larger octa¬ 
hedron (A3, panel (c)) are occupied with magnetic mo¬ 
ment 0.82 /i^ whereas the Ti states of smaller octahedron 
(B3, panel (d)) have negligible occupancy with magnetic 
moment 0.04 fiB- Compared with the magnetic moment 
of Ti in the middle of GTO (0.88 fis)^ we can consider 
the charge ordering as nominally one and zero electron 
for larger and smaller octahedron, respectively. The oc¬ 


cupied Ti atoms have ferro-orbital (in this case yz) and 
ferromagnetic order but with spin direction opposite to 
that in the GTO region. As pointed out in the previous 
section, the change in the orbital character from xz to 
yz along the A sublattice (panel (b) and (c)) is related 
to the change in the bond disproportionation between 
GTO region and interface octahedra. This in turn en¬ 
hances the localization of interface Ti states due to the 
small transfer integral between different t 2 g orbitals. We 
note that the occupied states at the larger interface oc¬ 
tahedron (panel (c)) have larger binding energy than the 
occupied states of GTO region (panel (a-b)) although the 
interface octahedral volume is slightly smaller than that 
of GTO region. Away from the interface the orbital and 
charge ordering disappears and the unoccupied t 2 g bands 
are shifted closer to the Fermi energy. The PDOS for 
other (GTO)^(STO)n interfaces with n > 1 are similar. 


V. ELECTRONIC STRUCTURES WITH A 
SINGLE SrO LAYER 


In this section, we present the electronic structure of 
(GTO) 5 (STO) 1 superlattices. In these superlattices the 
critical U for structural transition is significantly lower 
than in superlattices with multiple SrO layers. Unlike 
the previous cases with n > 1 , we will show that charge 
ordering is not a necessary condition for obtaining a in¬ 
sulating ground state and will show that the insulating 
phase persists to smaller U value due to the lack of low 
lying unoccupied states. 

The major difference between the electronic structures 
for single and multiple SrO layers is the absence, in the 
single-layer case, of the states coming from Ti atoms in 
between SrO layers. This significantly decreases the crit¬ 
ical of U for metal-insulator transition. Fig. 11 presents 
the DOS and PDOS of (GTO )5 (STO)i superlattices with 
U = 3 eV. From the total density of states we can see 
that the system is insulating with a band gap about 50 
meV. As seen in the PDOS, the Ti atom surrounded by 
GdO layers (panel (b)) has a density of states similar to 
that of bulk GTO. In our results with U = 3 eV there ex¬ 
ists both charge and orbital ordering but the difference in 
the magnetic moments of large (0.75 ys) and small octa¬ 
hedra (0.14 jub) is less than observed for U = 4 eV. This 
indicates that in-plane charge ordering is not crucial to 
the insulating phase. This conclusion is supported by the 
following argument. With two Ti 02 layers, yz (xz) or¬ 
bitals in each layer form bonding and antibonding states 
with one dimensional dispersion in y {x) direction. Due 
to this one-dimensional spectrum, if the electrons are fer- 
romagnetically ordered, the band is half-filled with strong 
Fermi surface nesting. Given that the inter-plane hop¬ 
ping is not significantly smaller than intra-plane hopping 
and xy band is higher in energy, a small orbital ordering 
that breaks the translation symmetry in both x and y di¬ 
rection and can open a gap at the Fermi energy. Thus as 
long as the octahedral distortion and electron correlation 
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FIG. 8: (color online) Projected density of states (PDOS) of Ti-t 2 ^ orbitals. The projection is done with t 2 g orbital defined 
with pseudo cubic axes that deviate with octahedral coordinates but the change is not significant. Blue solid lines, red dashed 
lines, and green dotted lines represent dxy^dyz, and dxz orbitals, respectively, (a) PDOS of GTO with U = 3 eV. (b-f) PDOS 
of (GT0)4(ST0)4 with = 3 eV for A-sublattice Ti atoms - numbers indicate atoms defined in Fig. 2. 



Ti02 plane 


FIG. 9: (color online) Orbit ally resolved and total occupancy 
of Ti-t 2 ^ states computed for = 3 eV (metallic regime) 
of (a) (GT 0 ) 4 (ST 0)4 and (b) (GTO) 4 (STO )6 super lattices. 
The gray region represents the Ti atoms between two GdO 
planes. 


keep the xy band unoccupied, the system becomes insu¬ 
lating even at weak correlation. We note that adding 
an additional STO layer will make the system metallic 
with U = 3 eV since isolating a lowest t 2 g band is much 


harder because the low lying unoccupied band from the 
added STO layer will cross the Fermi energy. In other 
words, lack of low lying unoccupied states and separa¬ 
tion of bonding and antibonding bands makes it possi¬ 
ble for n = 1 super lattices insulating for lower U value. 
Our results are consistent with the previously reported 
Hartree-Fock calculation by Chen et but for larger 
U we find sublattice charge ordering in additional to the 
orbital ordering. 


VI. PHASE DIGRAM FOR STRUCTURAL AND 
ELECTRONIC PHASES TRANSITIONS 


In the section, we present a phase diagram summa¬ 
rizing our computed results for structural and electronic 
phases for (GTO)^(STO)n superlattices. Fig. 12 shows 
the electronic and structural phase diagram in terms of 
critical U values and inverse of the STO layer thickness 
n. We can see that for n > 1, the MIT is accompanied 
by the structural transition which is crucial for the insu¬ 
lating phase. There is a narrow region of metallic phase 
before the structural transition with in-plane charge or¬ 
dering for n > 1 superlattices but the density of state 
is small at the Fermi energy. We identify this as a bad 
metal. 

Superlattices with a single SrO layer show qualitatively 
different phase behavior. Unlike the n > 1 case, there is 
no sharp change in the electronic structure across the 
structural transition from CO to CO II phase and the 
MIT occurs maintaining a weak layered charge order. 
This supports the idea^"^ that charge ordering is not nec¬ 
essary for n = 1 superlattices. Moreover the value of the 
critical U for the MIT of the n = 1 superlattice is signif¬ 
icantly smaller than n > 1 superlattices which also sup¬ 
ports the idea that weak correlation can open the band 
gap at the Fermi energy. Experimentally there is a MIT 
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FIG. 10: (color online) (a-f) PDOS of Ti-t 2 ^ orbitals for (GTO) 4 (STO )4 superlattice with eV following the atomic index 
defined in Fig. 2. The label at the bottom left of each panel represents Ti atoms where the letter A and B denote sublattice 
index and numbers indicate atoms defined in Fig. 2. Blue solid lines, red dashed lines, and green dotted lines represent dxy^ dyz, 
and dxz orbitals, respectively. 



E(eV) 

FIG. 11: (color online) Density of states of (GTO) 5 (STO)i 
superlattices with U = 3 eV. (a) DOS near the Fermi level, (b- 
d) Projected density of states of Ti-t 2 ^ orbitals for Ti positions 
defined in Fig. 4. Blue solid lines, red dashed lines, and green 
dotted lines represent dxy^dyz, and dxz orbitals, respectively. 



CO (Ins.) 
CO II (Ins.) 
Bad Metal 
Metal (I) 
Metal (II) 


FIG. 12: (color online) Electronic and structural phase dia¬ 
gram for (GTO)m(STO)n super lattices. We define in-plane 
and layered charge ordered insulating phase with orbital order 
as GO and GO II in gray and dark gray color, respectively. 
The structural GO phase is denoted as hatched lines so that 
the structural transition from GO phase to NGO (or GO II for 
n = 1) is expressed as the boundary of the shaded region. For 
n > 1, the narrow metallic region with charge ordering hav¬ 
ing very small density of state at the Fermi level is denoted 
as bad metal phase in light blue color. The ferromagnetic 
metal without charge ordering is denoted as metal phase in 
dark blue color and metallic phase with weak-layered charge 
ordering as metal II in purple. For the n = 1,3 superlattices 
phase boundaries are obtained with m = 2, the n — 2 super¬ 
lattice with m = 2, and n = 4, 6 super lattices with m = 4. 
The use of different m values is related to a small increase in 
the MIT phase boundary of n = 2 case. 


in between n = 3 and n = 2.^^ This is different from the 
prediction from the GGA+U phase diagram where the 
thickness dependent MIT occurs as the number of SrO 
layer changes from two to one. 


VII. BROKEN INVERSION SYMMETRY AND 
FERROELECTRICITY 

As can be seen from the Ti displacements in Fig. 2(d) 
or Fig. 3(c-d), in appropriate circumstances the superlat¬ 
tices we consider can develop non-centrosymmetric dis- 
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tortions leading to in-plane ferroelectric polarization. In 
bulk GTO, the Gd ions move substantially off center 
0.5 A) relative to the ideal perovskite position due 
to the the a~a~c^ octahedral rotation but the displace¬ 
ments alternate from layer to layer in the (001) direction 
so the effects cancel and there is no net polarization, and 
in particular the Ti atoms remain at center of the sur¬ 
rounding oxygen octahedron. Replacing a plane of Gd 
by Sr breaks the translation symmetry so that the polar¬ 
ization no longer cancels as previously suggested.The 
symmetry breaking also leads to an off-centering of the 
Ti atoms, suggesting an additional ferroelectric compo¬ 
nent. This shift of Ti atoms suggests that there can be 
ferroelectric moment as in PbTiOs. 

Fig. 13 shows atomic displacements of Ti and Gd 
atoms relative to the center of surrounding oxygen 
atoms for insulating (GTO) 4 (STO )4 (1/ = 4 eV) and 
(GTO) 5 (STO)i {U = 3 eV) superlattices. In both super¬ 
lattices, the shift is significant for near interface atoms. 
For out-of-plane direction, only interface Ti atoms are 
significantly displaced but the net moment is zero by the 
cancelation from the opposite contributions from two in¬ 
terfaces. For in-plane shifts, interface Ti atoms have the 
same displacement direction as that of Gd atoms next 
to the interface suggesting that the in-plane shift of Ti 
atoms is influenced by the movement of Gd atoms from 
the a~a~c^ octahedral rotation. The magnitude of the 
in-plane shift is about 0.07 A for Ti and 0.5 A for Gd with 
(GT0)4(ST0)4 superlattice and 0.03 A for Ti and 0.5 A 
for Gd atom with (GTO) 5 (STO)i superlattice. For com¬ 
parison, in PbTiOs the Ti and Pb atoms shift against the 
oxygen atoms about 0.26 A and 0.42 A, respectively.^^ 
Therefore we can see that the ionic shift is dominated 
by rare-earth atomic displacements. The amplitude of 
rare earth ion motion in our superlattices is compara¬ 
ble to that found in the canonical ferroelectric PbTiOs, 
but the amplitude of Ti atom motion in the superlattice 
is a about a factor of 5 smaller than the corresponding 
motion in PbTiOs- There are also differences in ferro¬ 
electric polarization between the superlattices with even 
and odd number of STO layers. When there is an even 
number of STO layers, the net ferroelectric moment is 
zero due to the opposite polarization direction of two in¬ 
terfaces. Whereas with an odd number of STO layers, 
there is an in-phase shift from interface Ti atoms giving 
a small upward polarization and there is relatively larger 
upward polarization from the displacements of rare earth 
atoms. As seen in Fig. 13 (d) the uncompensated dis¬ 
placement of Gd atom located between the sixth and the 
first Ti02 layer is larger than the displacement of Sr atom 
in between third and forth Ti02 layers. Combined with 
the difference in nominal ionic charge between Gd (3+) 
and Sr (2+), there is a net upward polarization with an 
additional contribution from a small displacement of Ti 
atoms in layer 3 and 4 (Fig. 13 (c)). Therefore we ex¬ 
pect that superlattices with an odd number of SrO layer, 
or more generally superlattices having two near interface 
Gd atoms displaced in the same direction may have a 




1 2 3 4 5 6 

Ti02 Layer 


FIG. 13: (color online) The displacements of Ti and rare earth 
atoms (Gd, Sr) for insulating superlattices, (a-b) The dis¬ 
placements in (GT 0 ) 4 (ST 0)4 superlattice with 1/ = 4 eV. 
(c-d) The displacements in (GTO )5 (STO) i superlattice with 
U — 3 eV. The displacements of Ti and rare earth atoms 
are defined as the distance from the center of mass of nearest 
neighbor in-plane oxygen atoms. The shaded regions repre¬ 
sents the bulk-like GTO region. The difference in the dis¬ 
placements between A and B sublattice is small so averaged 
values are shown. 


net ferroelectric moment. In the calculations performed 
so far, the polarization lies in the xy plane, but ‘tricolor’ 
superlattices with Gd-Sr-X may also have a 2 :-direction 
polarization. 
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VIII. BAND ALIGNMENT AND 
METAL-INSULATOR TRANSITION 

In the polar catastrophe scenario of LAO/STO in¬ 
terfaces, a MIT occurs as the number of LAO layers 
increases.The MIT is caused by a change in the band 
alignment of the STO conduction band driven by an 
increase in the electrostatic potential proportional to 
the LAO thickness.Similarly we can think of the 
metal insulator transitions of n > 1 superlattices in terms 
of band alignments between the filled lower Hubbard 
band of GTO and the low-lying conduction band of STO. 
Fig. 14 shows the PDOS of oxygen p and Ti d orbitals in 
the bulk-like GTO region, at the interface, and bulk-like 
STO region. For the insulating interface, we can see that 
the lower Hubbard band of GTO in the bulk-like region 
lies slightly lower than the band bottom of conduction 
band in the bulk-like STO region. On the other hand, for 
metallic interfaces, the conduction band of the bulk-like 
STO region overlaps the occupied lower Hubbard band 
of the bulk-like GTO region. The band overlap is deter¬ 
mined by four factors: the band gap in the GTO region 
Ag, Apd defined as the energy difference between oxy¬ 
gen p bands and the middle of upper and lower Hubbard 
bands, Asto denoting the band gap of STO, and Epoi 
defined as the energy shift in the top of oxygen p bands 
between GTO and STO region. We note that the energy 
shift Epoi arises from the dipole formed between GTO 
and STO. In the insulating phase, the energy shift is 
Ed/2 27re‘^n2Dd/e where £ is electric field between in¬ 
terface Ti02 and adjacent GdO layer, d is Ti-Ti distance, 
e is a dielectric constant, and n 2 D is the sheet charge den¬ 
sity of half electron per in-plane unit cell. In the metallic 
phase, electrons are delocalized over several layers of the 
STO region so the electric charge is farther from the inter¬ 
face and a larger dipole moment is expected. The energy 
shift is dz£{z) where L the distance from the 

GdO layer next to the interface to the Ti02 plane at the 
center of STO. The estimated energy shift for insulat¬ 
ing superlattice with four STO layers is about 155 meV 
by substituting the sheet charge density of 0.5 electron 
per in-plane unit cell and dielectric constant 75^^ and 
388 meV for metallic superlattice by assuming uniformly 
distributed electron gas (as in Fig. 9). These values are 
consistent with the energy shift in Fig. 14 which is about 
200 meV in the insulating superlattice and 700 meV for 
metallic superlattice. The factor of two differences may 
come from the uncertainty in the dielectric constant. 

Thus, we can write the condition for obtaining an in¬ 
sulating superlattice as follows: 

Epoi + ^STo{n) > Apd — Ag{U)/2 . 

Given that Apd and Epoi is not sensitive to the value 
of U and n, we can see that the critical U for MIT in¬ 
creases as the thickness of STO increases since Asto 
decreases as the number of STO layers n increases due 
to the confinement effect,while Ag is proportional to 
U. Although we neglect the energy cost of octahedral 



FIG. 14: (color online) PDOS of Ti d (blue solid line) and 
oxygen p (dotted red line) orbitals of (GTO) 4 (STO )4 super¬ 
lattice. (a) PDOS in the middle of GTO with 1/ = 3 eV and 
(b) 1/ = 4 eV. (c) PDOS of interface Ti02 layer for 1/ = 3 eV 
and (d) U=4 eV. (e) PDOS in the middle of STO for U=3 
eV and (f) U=4 eV. The solid black line is the Fermi energy 
and the dashed black lines in the panel (a) and (b) represent 
the band edge of oxygen p orbital in the middle of GTO. The 
energy difference between the dashed black line and dotted 
black line (center of gap between lower and upper Hubbard 
bands) is defined as Apd- In panel (e) and (f) the Epoi is 
defined as the energy difference of oxygen band edge of GTO 
and STO region. 


distortion and energy gain from localizing electrons, this 
relation gives a reasonable estimate of the ground state of 
heterostructures between band and Mott insulator with 
bulk properties the materials. 
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IX. SUMMARY 

We investigated the metal-insulator transition of 
GTO/STO superlatives using first-principles GGA+U 
method. Two different mechanisms for the insulating 
phase are identified. We find that charge and orbital or¬ 
dering accompanied by difference in the volume of inter¬ 
face octahedra (“two sublattice charge order”) are neces¬ 
sary to obtain the insulating phase for superlattices with 
n > 1. On the other hand, in superlattices with n = 1 
the insulating gap emerges via a combination of orbital 
ordering (leading to a quasi-one dimensional in-plane dis¬ 
persion) combined with a bonding-andtibonding splitting 
arising from coupling across the SrO layer, as previously 
found by Ghen et We find that the critical U needed 
to drive the MIT for n = 1 superlattices is significantly 
smaller than n > 1 since small orbital ordering can open a 
gap in the bonding band due to the instability in the one 
dimensional dispersion of yzjxz bands. A local inversion 


symmetry breaking around Ti atoms is observed and it is 
shown that ferroelectric polarization is possible with odd 
number of STO layers. We present the phase diagram 
for general (GTO)m/(STO)n and within GGA+U show¬ 
ing discrepancy between our results and the transport 
measurement of Ref. 37, which shows a metal insula¬ 
tor transition already at n = 2. The difference cannot 
be resolved by fine-tuning U implying the importance of 
additional physics not included in the DFT+U approxi¬ 
mation. 
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